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Abstract: We revisit the familiar construction of one-loop scattering amplitudes 
via generalized unitarity in light of the recently understood properties of loop inte- 
grands prior to their integration. We show how in any four-dimensional quantum 
field theory, the integrand-level factorization of infrared divergences leads to twice 
as many constraints on integral coefficients than are visible from the integrated ex- 
pressions. In the case of planar, maximally supersymmetric Yang-Mills amplitudes, 
we demonstrate that these constraints are both sufficient and necessary to imply the 
finiteness and dual-conformal invariance of the ratios of scattering amplitudes. We 
present a novel regularization of the scalar box integrals which makes dual-conformal 
invariance of finite observables manifest term by term, and describe how this proce- 
dure can be generalized to higher loop-orders. Finally, we describe how the familiar 
scalar boxes at one-loop can be upgraded to 'chiral boxes' resulting in a manifestly 
infrared-factorized, box-like expansion for all one-loop integrands in planar, M = 4 
super Yang-Mills. Accompanying this note is a Mathematica package which im- 
plements our results, and allows for the efficient numerical evaluation of any one-loop 
amplitude or ratio function. 
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1. Introduction 

One-loop amplitudes have been extensively studied in recent decades, leading to 
many important insights and discoveries about the structure of scattering ampli- 
tudes, and frequently serve as an important source of theoretical 'data' with which 
to test new ideas [1-11]. A powerful approach to computing loop amplitudes in 
any quantum field theory is the unitarity-based method, in which the amplitude is 
expanded into a basis of standardized scalar Feynman integrals (regulated if neces- 
sary) with coefficients fixed by on-shell scattering processes. Although very familiar 
and reasonably well understood, the way this approach has been realized in terms 
of existing technology does not make manifest several recently discovered aspects of 
loop-amplitudes — especially for the particularly rich case of scattering amplitudes in 
planar, N = 4 super Yang-Mills (SYM). 

The two principle shortcomings about the way generalized unitarity has been 
realized in terms of existing tools (at least for N = 4 SYM) are: (1) that it fails to 
reflect the rich symmetries observed in loop amplitudes prior to integration; and (2) 
that even those symmetries which survive integration — such as the dual-conformal 
invariance (DCI) of the ratios of scattering amplitudes (see e.g. [12]) — are severely 
obfuscated in all existing regularization schema for infrared-divergent contributions. 
Because of this, manifestly-DCI expressions for ratio functions are known only in a 
few exceptionally simple cases (see e.g. [12-16]). In this note, we revisit this story and 
fully address both shortcomings, providing manifestly-DCI expressions for all one- 
loop ratio functions in M = 4 SYM, and describing how the familiar box expansion 
can be upgraded to a chiral box expansion which matches all one-loop integrands. 

This paper is organized as follows. In section 2 we review how generalized unitar- 
ity can be used to reproduce integrated amplitudes in J\f—4 SYM, and in section 2.2 
we (heuristically) derive 'D CP-regularized expressions for all scalar box integrals, 
which are given in Table 1. In section 2.3 we summarize the computation of scalar 
box coefficients using momentum-twistor variables, and write explicit formulae for 
all one-loop box coefficients in Table 3. 

In section 3 we explore the general features of the 'DCI'-regularization proposal. 
In section 3.1 we show that this proposal correctly reproduces all finite observables of 
any planar theory, thereby justifying its description as a 'regularization scheme'. In 
section 3.2 we describe how this scheme can be extended beyond one-loop, and com- 
pare it with existing approaches. The 'DCI'-regulator is closely related to (and moti- 
vated by) the way that infrared divergences arise at the level of the loop-integrand. In 
section 3.3, we describe how the IR-divergences of loop amplitudes appear in terms 
of the 'DCI'-regularization scheme, and in section 3.5 we show how generalized uni- 
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tarity realized at the integrand-level can be used to generate more powerful identities 
than would be possible after integration. In section 3.6 we illustrate the how these 
features persist beyond the planar-level. 

In section 4 we return to the familiar box expansion, and describe how it can be 
made 'chiral' in a natural way, allowing us to match the full, chiral integrand of any 
one-loop scattering amplitude in J\[=&. 

For the purposes of concreteness and completeness, we review the basic kine- 
matical variables — momentum twistors — used in most of this paper in Appendix A; 
and in Appendix B, we use these variables to give a closed- form specialization of the 
BCFW recursion relations for all one-loop integrands in A/" =4 SYM. 

And finally, we have implemented the results described in this paper in a Math- 
ematica package called 'loop_amplitudes' which is documented in Appendix C. 
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2. Revisiting Generalized Unitarity at One-Loop 



A major triumph of the unitarity-based approach to quantum field theory was the 
discovery that any one-loop amplitude can be written as a linear combination of stan- 
dardized, scalar integrals with coefficients expressed as on-shell diagrams, historically 
known as 'leading singularities' (for a comprehensive review, see [11]). Because of 
the good UV-behavior of A^ = 4 SYM, only box integrals contribute — those involving 
four loop- momentum propagators — giving rise to the familiar 'box expansion': 




c C c C 

(To be clear, throughout this paper we will refer to (Jd 4 £)An^' 1 as the (integrated) 
one-loop N fc MHV amplitude, expressed in units of g 2 N c / (16n 2 ).) 

The objects appearing in (2.1) will each be described in detail below. But let 
us briefly remark on the motivation underlying the box expansion. Loop amplitudes 
are obtained by integrating the loop integrand over a four- dimensional contour of 
real loop momenta, I G 1R 3 ' 1 . If this integrand were obtained from the Feynman 
expansion, for example, it would include many propagators for the internal, 'loop' 
particles. A co-dimension one residue of this integral 'enclosing' one propagator — a 
'single-cut' — would correspond to putting one internal particle on-shell. Because the 
loop integral is four-dimensional, the highest-degree residues would be co-dimension 
four; these are the so-called leading singularities or 'quad-cuts' of the integrand 1 and 
are computed in terms of on-shell diagrams of the form shown in (2.2). 

The scalar box integrals I a ,b,c,d are simply those involving precisely four Feynman 
propagators of a scalar field theory — normalized to have co-dimension four residues 
of unit magnitude. As such, we should be able to represent any integrated loop 
amplitude in M = 4 by dressing each box integral with the actual co-dimension four 
residues 'enclosing' the corresponding propagators as for the full loop integrand. A 

1 In a general quantum field theory, the integrand may also have contributions where only co- 
dimension three (triangles) or co-dimension two (bubbles) residues exist; it is a non-trivial fact that 
TV = 4 does not need contributions from triangles or bubbles (closely related to its UV-finiteness) . 
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slight subtlety is that the residues of scalar boxes come in parity conjugate pairs, so in 
order to agree with the complete integrand the scalar boxes should be supplemented 
by parity-odd integrals, [8]. Since these integrate to zero, they are often ignored. 
For the purpose of computing the integrated amplitude in this section — as opposed 
to the integrand — we will also ignore them here. This mismatch will be addressed in 
greater detail in section 4, where we show how to upgrade the box expansion (2.1) 
in a way which allows us to match the full amplitude prior to integration. 

2.1 Scalar Box Integrals and their Divergences 

Let us start our analysis with the generic, 'four-mass', scalar box integral [17,18]- 
one for which all four corners are 'massive' (involving at least two massless momenta): 
A a 



la 



b,c,d 



d 

D 



4¥ 



B _ [« B -(a,c)(b,d)A 



b J (£,a)(£,b)(£,c)(£,d) 



Jd^£Xa, b ,cd{£), (2.3) 



where 



A^(l-u-„)>-4™, (2.4) 

(a,c)(b,d) (a,c){b,d) 

with p a = x a+ i-x a [19] and where l/(£,a) denotes the standard propagator, 

(a,b) = (x a -x b ) 2 = (p a + Pa+i+ ■ ■ ■ +Pb-i) 2 and (£, a) = (£-x a ) 2 . (2.5) 

When all the corners are massive, this integral is a transcendentality-two, com- 
pletely finite, symmetric function of the dual-conformally invariant cross-ratios (u, v): 



-I a ,b,cA u i v ) = Li 2(«) + Li 2 (/3) - Li 2 (l) + - log(n) log(v) - log(a) log(/3), 



(2.6) 



with 

a = -(l-u + v + A) and (3 = -(1 + u - v + A). (2.7) 

Notice that this form is manifestly symmetric under the exchange u -h- v as this 
exchange results only in a -H- /3, under which (2.6) is obviously symmetric. The 
equivalence of (2.6) to existing formulae in the literature is easily verified 2 . 

It is easy to see that this integral becomes divergent when any corner becomes 
massless — for example, identifying legs p a and ps results in: 

2 For notational simplicity, we have left implicit a factor ir 2 in the definition of the measure of 
the integrand (2.3) — the integrated expression, (2.6), is of course the standard one; we also omit 
factors of i associated with the Wick rotation: "<i 4 £" denotes the Euclidean integration measure. 
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c C 

This causes the cross ratio u to vanish, introducing a logarithmic singularity from 
the term — -= log(w) log(t>) in (2.6). It is worth noting that (2.6) simplifies consid- 
erably when u is taken to be parametrically small, u — > 0(e): A — > (l—v) + 0(e), 
a-)-l + 0(e), and (1 -v ) + C(e), resulting in 

-iaAcdM) = Li 2 (l-^ + -logHlog(^ + C>(e). (2.9) 

M-s>C(e) Z 

This divergence can be regulated in a number of ways, including dimensional regular- 
ization (see e.g. [18] for standard formulae). Another canonical way to regulate such 
divergences uses the Higgs mechanism, [20]. In the simplest setup, this is used to 
give masses only to internal propagators, and leads to the mass-regularized formulae 
found in e.g. [21,22]. However, because such regularization schema are predicated 
on a dimensionful parameter, the regulated formulae that result severely break dual- 
conformal invariance, obscuring the ultimate invariance of even absolutely convergent 
(and hence DCI) combinations of box integrals. (A complete basis involving only 
manifestly finite integrals for all convergent one- loop integrals was given in ref. [16].) 

We are therefore motivated to find some way to make all the singular limits of 
the general integral (2.6) as dual-conformally invariant as possible, regulating the 
divergence caused by u — >■ in a way which depends only on some dimensionless 
parameter, denoted e, and dual-conformal cross-ratios. Such a regularization scheme 
is described in the next subsection. 

2.2 The 'DCI' Regularization of Scalar Box Integrals 

We propose the following regulator for one-loop integrals: render all external legs 
slightly off-shell by displacing the coordinates according to, 

- / \ (a-2, a) /n-irA 

x a ^x a = x a + e(x a+1 -x a )- — ; (2.10) 

{a -2, a+1) 

this transformation can be understood graphically as follows: 
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We will prove that this regulator produces the correct result for all finite observ- 
ables and discuss its implications, generalizations, and extensions in greater detail 
in section 3; but first, let us understand its consequences for the one-loop integrals 
appearing in the box expansion (2.1). 

If the dimensionless parameter e is small (the only regime in which we will be 
interested) then the invariants (a, 6) which are already non-vanishing are modified by 
a negligible amount. However, when b = a+l as in (2.8), for example, the invariants 
will be regularized by (2.10): 

~ (o-l^+l) + w . th 6 = fl+i _ 

(a- 1,6+1) 

This implies that the integrated expressions for all box integrals will be given by 
limits of the four-mass-box (2.6); in particular, it is not necessary for us to add any 
'discontinuity functions' such as those discussed in e.g. [4]. 

Consider the leading, degenerate limit of the box function, the so-called 'three- 
mass' integral. This corresponds to taking 6 = a+l in (2.6), keeping all other external 
corners massive: (6, c), (c, d), (d, a) ^ 0. The cross-ratio v is unaffected by the regu- 
lator to 0(e), while u transforms according to 

(a,c)(b,d) (a-l,b+l)(a,c)(b,d) 
Using (2.9), the limit is easily seen to be given by 
^ a 



d 
D 



4/ 



-Li 2 (l - v)- 1 - log(u') log(u) - l - log(e) log(v) + 0{e). (2. 13) 



c C 

All further singular limits of (2.6) are similarly regulated according to the shifts 
(2.10). This results in 'DCI-regularized' forms for all scalar box integrals, which 
we have listed in detail in Table 1. Importantly, these integrals depend only on the 
dimensionless parameter e and dual-conformally invariant cross-ratios. 

One special case not given in Tablel is the so-called 'massless' scalar box — which 
is only relevant for n = 4. In this case, the shifts (2.10) are not strictly defined, but 
can be generalized in simple way so that u,v— >e 2 , resulting in the following, 'DCT- 
regulated massless box function: 

-/ 1)2 , 3i 4 = Li 2 (l) + 2 log(e) 2 + 0(e) (only for n = 4). (2.14) 
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Table 1: 'DCI'-Regulated Scalar Box Integrals 



A a 






6 C 


-I = Li 2 (a) + Li 2 (/3) - Li 2 (l) + - Iog(«) log(i;) - log(a) log(/3) 


(a,b)(c,d) (b,c)(a,d) , 
^ ^ , ' ^ , and = §(1 +u -« + A 
(a,c)(M) (M(M) A = ^/(l -u-t>) 2 -4uf 


A n 
D — ^ ^ b 


-F = Li 2 (l - t/) + i log(w') log(«) + i log(e) log(«) + 0(e) 


c C 


u I— > en' 


, _ (a - 1, 6)(o, 6+1) (c, d) (6, c)(a, ti) 


ti = a 


V I— >■ f 


(a-l,b+l)(a,c)(b,dy ' ~~ {a lC ){b,d) 


A a ^fi 




A 2 (l) + 1 log(uO log(^) + i log(e) log(«V) + ^ log(e) 2 + 0(e) 




-I £ = I 


UH-> Ell' 


/_ ( a -l^)(c,rf) ,_ (6,c+l)(a,d) 


d = a, G = o 




(a-l,c)(b,dy ' ~~ (b,d){a,c+l) 


A a 


-I e = Li 2 (l - v) + i log(u') log(^) + log(e) \og{v) + 0(e) 


c c 






u h+ e 2 it' 


,_ (a-l,b)(a,b+l)(c-l,d)(c,d+l) _ (b,c)(a,d) 


n = a, L) = c 




~ (a-l,b+l)(c-l,d+l)(a,c)(b,d)' 1 ~ {a,c)(b,d) 


A a 

"ff 


-I e = Li 2 (l) + 1 log(e) log(u') + log(e) 2 + 0(e) 








u I-+ e 2 it' 


, _ (a-l,b)(c,d+l) 


B = a, C = b, D = c 




(a-l,c)(b,d+l) 
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2.3 Scalar Box Coefficients: One-Loop Leading Singularities 

Factorization dictates that the residues of the loop integrand — called leading singu- 
larities — are simply the products of tree-amplitudes, summed-over all the internal 
particles which can be exchanged, and integrated over the on-shell phase space of 
each. As mentioned above, we represent such functions graphically as on-shell dia- 
grams of the form shown in (2.2). 

These are simply algebraic (super)functions of the external kinematical data — 
almost always rational, and at one-loop involving at most the solution to a quadratic 
equation. Leading singularities have of course been known to the literature for quite 
some time, and can be computed in many ways. A comprehensive summary of 
these objects — their classification, evaluation, and relations — was described recently 
in ref. [11]. The physical content of any on-shell diagram (after blowing up all tree- 
amplitudes at the vertices themselves into on-shell diagrams) is encoded by a permu- 
tation [11], and the permutations of the corners are simply 'glued' together to give 
the permutation of the 'one-loop' diagram: 




(2.15) 



Given the permutation labeling an on-shell diagram, it is trivial to construct an 
explicit formula for the corresponding on-shell function. This is done most directly 
in terms of an auxiliary Grassmannian integral as described in ref. [11]. (All the 
necessary tools involved in this story have been made available in a Mathematica 
package called 'positroids', which is documented in ref. [23].) We will not review 
these ideas here, but simply give the formulae that result. 

The most compact expressions for leading singularities are found when they are 
written in terms of the momentum-twistor variables introduced in ref. [24] since 
they simultaneously trivialize the two ubiquitous kinematical constraints — the on- 
shell condition and momentum conservation. Momentum-twistors are simply the 
twistor-variables [25] associated to the region-momentum coordinates x a defined in 
section 2. (A brief introduction to momentum-twistor variables and an explanation 
of the notation used throughout this section is given in Appendix A.) 

Assuming a modicum of familiarity with momentum-twistors, let us now describe 
the form that leading singularities take. We start with the most general 
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leading singularity involving four massive corners. It turns out that this case is the 
only one we need to consider, as it will smoothly generate all the others in a very 
natural way. 

The most important data are the two solutions l\,£ 2 to the kinematical con- 
straints of putting all the internal lines on-shell, 



A a 




(2.16) 



Here, the lines (Aa), . . . , (Dd) in momentum-twistor space correspond to the region- 
momenta x a , . . . ,Xd used in section 2.1; and the lines £\ and £2 correspond to the 
two solutions to the problem of putting four propagators on-shell. These 'quad-cuts' 
are the solution to a simple geometry problem in momentum-twistor space (viewed 
projectively as P 3 ): £\ and £2 are the two lines which simultaneously intersect the 
four generic lines (Aa), . . . , (Dd): 




(2.17) 



(The fact that there are two solutions to the problem of putting four-propagators 
on-shell is a classic result of the Schubert calculus — and continues to hold even when 
the four lines are non-generic; see ref. [16] for an exposition of these ideas.) 

We will give explicit formulae for the twistors cti, . . . , 5\ and a 2 , ■ ■ ■ ,5 2 corre- 
sponding to the two lines t\ and £ 2 intersecting the lines (Aa), . . . , (Dd), respectively, 
in Table 2; but for now let us take it for granted that they are known. Given these 
twistors, it is easy to write the leading singularity for each particular quad-cut £i y 2 
in terms of momentum-twistors. In momentum-twistor space, we are dealing with 
MHV-stripped amplitudes (so (N fc=0 )MHV tree-amplitudes are simply the identity), 
and polarization sums become simple multiplication of the corresponding MHV- 
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stripped amplitudes. This allows us to 'peel-off' the tree-amplitudes at each corner 
from a standard on-shell graph involving only MHV amplitudes, [13]: 

A a A a 



( A a {ai,a,. 
xAb(Pi,b,. 
xA c (ji,c, . 



c C c C 

where the on-shell graph on the right is the N fc=2 MHV 'four-mass' function 3 [11], 
A a 



c C 

While simply replacing (ai, . . . , 5{) — > (« 2 , . . . ,62) in the formula above would give 
(minus) the other leading singularity, we will find it advantageous to use an alternate 
form of the four-mass function involving the £2 solution: 

A a 





[a 2 BbCc}[j 2 DdAa] 1 



(a 2 cDd){^2 & Bb) 
(a 2 cBb)(^ 2 aDd) 



(2.19) 



In Table 2, we give the particular solutions to the quad-cuts represented graph- 
ically in (2.17) — where A is defined as in equation (2.4). (The notation used here is 
fully defined in Appendix A.) The motivation for using two separate formulae for the 
four-mass leading singularities is that they separately evolve smoothly to all other 



3 We note that this expression for the four-mass function appears to differ from that given in 
ref. [26] which we suspect to be incomplete. 
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cases. This is made possible by the fact that the multiplicative factors appearing in 
Table 3, which encode the shifts of the quad-cuts from each 'corner' of the box (see 
equation (2.17)), are all smooth and non-singular in limits where some of the legs 
are identified. (Notice that (pi j2 — ^1 i n all limits where a pair of legs are identified.) 

As promised, the formulae given above for the four-mass leading singularities 
smoothly generate all one-loop leading singularities, which for the purposes of com- 
pleteness and reference have been written explicitly in Table 3. The formulae given 
in Table 3 have been organized in order to highlight how each case descends smoothly 
from the general cases given above for the four-mass functions. 

The complete box coefficient for a given topology is given by the sum of the two 
corresponding on-shell diagrams, f a ,b,c,d = fabcd Jr fabcdi which involve the quad-cuts 
l\ and £ 2 , respectively: 



A a A a 




c C c C 



Here, each graph represents a sum over all such graphs with the same topology — 
involving any four amplitudes An^\ . . . ,Au4 such that k a + kb+ k c + kd = k-2, and 
for which n a + nb+n c + n<i = n+8, with <k < n-4 for each (except when n = 3, for 
which A^f 1 ^ is allowed). 



, , x ~,„, , {aBb(cC)n(DdA)) + {ABb(cC)f)(Dda)) + {aAcC){BbDd)A 



Pi=(Bb)r\(aA6i)=z B +Zb 



2(Bb {cC)O(DdA)A) 
(B aA (Dd)f](cC b)) + (b aA (Dd)f](cC B)) + (aA cC) (Bb Dd)A 



2(aA (Dd)f)(cCb) b) 
7i = (cC)n(-Ddai) and 8\ = (Dd)f){cC fa) 



^rv^n W Z (AdD(Cc){)(bB a)) + (adD(Cc){)(bBA)) + (AaCc)(bBdD)A 



p 2 =(bB)r\(Cc5 2 ) =z b +z B 



2{dD(Cc)f]{bBa)a) 
(b Cc (dD)f](Aa B)) + (B Cc (dD)f](Aa b)) + (Aa Cc) (bB dD)A 



2(Cc (dD)f)(Aa B) B) 
72 = (Cc)C\(bB a 2 ) and 5 2 = (dD)f)(Aa (3 2 ) 



Table 2: Explicit solutions £1,^2 to the Schubert problem involving four generic lines. 
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Table 3: All One-Loop Leading Singularities in Momentum-Twistor Space 




</>i[<Si AaBb]\PiCcDd]x. 
A a {on, ■ ■ . ,/3i)A(/3i, ■ ■ ■ ,71) 
A c (ji, . . .,8i)Ad(5i, ...,ai) 

a 1 ^(aA)f](Bb- fl ) \ 
/3 1 =(Bb)r\(aA6 1 ) I 
7i=(cC)n( J Drfa 1 ) | 
S 1 =(Dd)r\(cC/3 1 )) 



(p 2 [a 2 BbC c] [-f 2 DdAa]x 
A a {a 2 , . . • ,/3 2 )A(/?2, • ■ ■ ,72) 
A(72, ■ ■ ■,S 2 )Ad(.S 2 , ...,a 2 ) 
a 2 =(Aa)f]{dDj 2 ) 
2 =(bB)f]{CcS 2 ) 
l2 = {Cc)r\(bB a 2 ) 
S 2 =(dD)r\(Aap 2 ) 





B [fcCcDdjxAbiPu---,!!) 
-4 c (7i, • ■ ■ ,5i)A d (Si, . . . ,ai) 



[a 2 BbCc][j 2 DdAa]x 
A b (,B 2 ,... }l2 ) 
A c (j 2 , ■ ■ -,S 2 )Ad(S 2 , ...,a 2 ) 



ai—a 

H=(cC)f]{Dda) 
S 1 =(Dd)r\(.cCB) 



a 2 =(Aa)r\(dD l2) 
, P 2 ={bB)f\{Cc5 2 ) 
l2 =(Cc)pl(bBA) 
5 2 =(dD)f)(Aab) 




c C 




B [frCcDd]* 
A c (ji, ■ ■ ■ ,Si)Ad(Si, ■ ■ . 



a\—a 
0i= B 



yi =(cC)r\{Dda) 
C=h 1 S^iDd^icCB) 



[y 2 DdAa]x 
A(72, ■ ■ -,S 2 )Ad(.S 2 , ...,a 2 ) 

a 2 =(Aa)f)(dDC) 
P 2 =b 

12=C 

S 2 =(dD)r\(Aab) 




C=b 




A b (f3 1 ,...,ji)A d (8 1 ,. . .,ai) 



[a 2 BbCc][j 2 DdAa]x 
A b {(3 2 , . . .,~f 2 )A d (S 2 , ...,a 2 ) 



li- 



ly 



Pi=B 

7i=c 
6x=D 

A d (5 u . . . ,ai) 

ai—a 
A = B 
71= c 



a 2 =(Aa)n(dDC0 
; ^ 2 =(6B)n(Cai) 

72 =(cc)n(6s^) 

^=(rfg)f|04a6) 
[72-Ddyla] xA d (S 2 , ...,a 2 ) 



( a 2 =(Aa)f)(dDC) 
fo=b 
72= C 
{ 6 2 =(dD)C\{Aab) 



2< 




with 



-1 



l ~ {fiidCc)(5ibAa)) and ^ 2 



1 



(a 2 cDd){y 2 aBb) 



(a 2 cBb)(j 2 aDd) 

(Cyclically-related leading singularities are related by fl bcd = f bcda -) 



-1 
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3. Properties and Extensions of 'DCI' Regular izat ion 



Using the 'DCI'-regularized box integrals bcd and the box coefficients described in 
section 2, the scalar-box expansion becomes 

d^i A^' 1 = \fa,b,c,d + fa,b,c,d) ^a,b,c,d- (3-1) 

J \<a<b<c<d 

Let us now prove that this produces the correct result for all finite observables in 
planar Af = 4 SYM by showing that the expressions given for I e ahcd can be obtained 
from a very concrete and simple regularization procedure. 

3.1 The 'DCI' Regularization Scheme 

Given any four-dimensional integrand T(£), we define its 'DCI'-regulated integral by 
deforming the integrand in the following way: 

IdH X(£) = [d 4 £ 1(£)R(£) where R(£) = T\t^- (3-2) 
J J f = i( £ > a ) 

reg 

The regulating factor R(i) suppresses all IR-divergent regions (for all integrands), 
making the result manifestly IR- finite. Ultimately, this works because all infrared sin- 
gularities in a planar integral arise from predetermined integration regions — namely, 
the collinear regions f->as a +(l-a) x a +i, where the factor R{£) becomes 0(e). 

It is trivial to see that the regulated expression (3.2) produces the correct result 
for all finite observables in the limit of e— s-0: since the factor R(£) ~l + (9(e) every- 
where except in the isolated regions responsible for infrared divergences (where it is 
(9(e)), it can be ignored for any convergent integral. 

We should emphasize a distinction we are making between convergent integrals 
and so-called "finite" integrals. Tautologically, a "convergent" integral is one which 
can be evaluated without regularization; this requires that it have no collinearly- 
divergent regions. This notion of convergence precludes the combinations of divergent 
integrals which happen to be "finite" in some particular regularization scheme (but 
for which integration does require regularization). The convergence of an integral 
can be tested as follows: multiply the integrand by any two adjacent propagators, 
and verify that the product vanishes in the corresponding collinear region: 

lim [(£,a)(£,b)xX(£ = ax a + (1 -a)x b + 5£)} = for b = a+1. (3.3) 

As pointed out in ref. [16] the integrand for the ratio function is in fact convergent. A 
similar integrand-level test of (partial) convergence for the logarithm of the 4-particle 
amplitude has been shown sufficient to completely fix the integrand through seven- 
loops, [27,28], and has also been used to find amplitudes involving more particles [29]. 
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It remains for us to show that the regulated amplitude (3.2) is actually given by 
(3.1). To see this, consider the box expansion as an integrand-level statement: we 
can decompose any one-loop integrand into parity-even and parity-odd sectors: 

4P' 1 V) = Yl (/i 6> c, d + /a 2 b ,M)^,c d W+paxity-odd. (3.4) 

l<a<b<c<d 

(Recall that T a ,b,c,d{£) denotes the integrand of the scalar box I a ,b,c,d-) The scalar 
boxes form a complete basis for parity-even four-dimensional integrands (see e.g. 
[18,30,31]), and so the first term of (3.4) completely captures all parity-even contri- 
butions. The parity-odd contributions in (3.4) are often ignored because they vanish 
when integrated over the parity-even contour of M 3 ' 1 . Importantly, all parity-odd 
integrals are not merely vanishing upon integration, but are in fact convergent in 
the sense described above, [16]. (This is not too surprising since the requirement for 
an integrand to vanish in the limit (3.3) is itself parity-invariant.) And because all 
parity-odd integrands are convergent, the regulator R(£) is l + 0(e) everywhere, and 
can therefore be ignored. 

For the parity-even sector — that is, the scalar box expansion, (2.1) (but now 
understood at the integrand-level) — we only need to verify the following: 

^ jf t C) }t d = fa ~M (M)A ~ + 0(e) = I lbcd . (3.5) 
(£,a)(£,b)(£,c)(£,d) J (£,a){£,b)(£,c){£,d) 

This identity is proved by noting that all other factors in R(£) are approximately unity 
except in very small regions, where the unregulated box is not singular for lack of 
divergent propagators; with the explicit propagators regulated, the singular regions 
are all removed, resulting in precisely /^ C(i as described in section 2.1 and given 
in Table 1. This concludes our proof that (3.1), using the 'DCI'-regulated scalar 
box integrals, correctly reproduces all regulator-independent contributions to loop 
amplitudes, and therefore leads to correct formulae for all finite one-loop observables. 

3.2 Generalization of 'DCI'-Regularization to Higher Loop-Orders 

The failure of simple off-shell regularization beyond one-loop has been known for 
quite some time (see e.g. [19]). For example, for the two-loop 4-particle amplitude in 
J\f—4 SYM, a simple off-shell prescription fails even to give the correct coefficient for 
the double-logarithmic divergence! Because the 'DCI'-regulator (3.2) is somewhat 
similar to an off-shell regulator at one-loop (but one involving non-uniform masses), 
this may seem like a bad omen for extending it beyond one-loop. However, the most 
natural generalization of (3.2) beyond one-loop cannot be interpreted as an off-shell 
regulator — which is good news, indeed! Let us first describe the generalization of 
the 'DCI'-regulator to higher loop-orders, and then illustrate how it differs from an 
of-shell regulator in the case of 4-particles. 
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The integrand-level understanding of the 'DCI'-regulator, (3.2), provides an ob- 
vious generalization to higher loops-orders, 

]Jd%l(£ 1 ,... } £ L )= ]J{d% R(£i)) I(£ 1 ,...,£ L ). (3.6) 

We suspect that this regularization prescription will render all multi-loop integrands 
IR-finite; moreover, using the same arguments as in section 3.1, we expect that (3.6) 
will generate the correct result for all finite multi-loop observables. 

As one simple example of a regularized integral beyond one-loop — and as an 
illustration of the differences between (3.6) and off-shell regularization — consider the 
'DCI'-regularized 4-particle double-box integrand, [32]: 

d%d u (1,3)(2,4) (4,4)(4,2) 



(4,1)(4,2)(4,3)(4,4)(4,3)(^,4)(^,1) (4,4)(4,2) 

Because of the part of i?(£i)i?(£ 2 ) which survives, this is clearly not an off-shell 
version of the double-box! In particular, the integrand (3.7) includes regions (both 
collinear and soft-collinear) where the regularization-factor cannot be approximated 
by unity. This is good news, since this had to happen for the regulator to have any 
chance of working beyond one-loop. 

It would be interesting to compute this integral explicitly, and verify for example 
that scheme-independent quantities such as the two-loop double-logarithmic diver- 
gences (the so called cusp anomalous dimension) are correctly reproduced; we leave 
this, however, to future work. 

3.3 The Infrared-Divergences of One-Loop Amplitudes 

Let us now consider how the IR-divergences of scattering amplitudes are organized 
in the 'DCI'-regulated box expansion of (3.1). These divergences arise from the parts 
of Iabcd proportional to log(e) or log(e) 2 ; let us denote the combined coefficients of 
each of these divergences as follows: 

JdH A^' 1 = -F 2 log(e) 2 - Ft log(e) - F + 0(e). (3.8) 

It is easy to identify the coefficient F 2 from Table l: 4 the only integrals which 
include a factor of log(e) 2 are the so-called 'two-mass hard' and 'one-mass' boxes, 

4 We should note that the coefficient of log(e) 2 in the amplitude is twice the coefficient of log(/i 2 R ) 2 
of the Higgs regulator described in ref. [20] — equivalently, it is twice the coefficient of 1/(2(D — 4) 2 ) 
appearing in dimensional-regularization. This can be understood from the way (3.2) cuts-out the 
soft-collinear divergent regions of integration. 
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F *= E (/U+V + C.aw)^^' (3-9) 



a+2<d<n 
15 



which can be understood as the 'BCF' representation of tree-amplitudes discovered 
in ref. [33] (see also [34]). 

Interestingly, the coefficient Fx in (3.8) is also proportional to the tree- amplitude. 
(This follows from the ideas discussed in section 3.5.) Because the coefficient of 
log(e) must be of transcendentality-one, it should involve the logarithm of some 
dual-conformal cross ratio, which we denote Q n ; ultimately, Fx is found to be simply, 

F 1 = 2^ log(^) with ^^TTr^l- ( 3 - 10 ) 

±± (a, a +3 1 

An important object that arises in the discussion of N fc MHV amplitudes A^' e 
is the so-called ratio function, denoted 7^1 , which at one-loop is given by, 

ni k)A = A m_ A (k),o xA m_ (311) 

Because the divergences Fx and F2 are always proportional to the tree-amplitude, 
An^'° — and because in momentum-twistor space An^'° is simply the identity — we see 
that IZ^' 1 is finite in the limit e— s-0. Therefore, any one-loop ratio function can be 
computed more simply as, 



f^f V (k),l_ v ff 1 + f 2 _ jw.o x f (MHV) V fin 

UZ IX, n — \Ja,b,c,d ' J a,b,c,d -^n x J a,b,c,d J 1 a,b,c,d^ 



l<a<b<c<d 



(3.12) 



where f^bcd are ^ ne MHV one-loop box coefficients, and 1^1 cd are the e-independent 
("finite") parts of the e- regulated scalar box integrals I^bcd as listed in Table 1. 
Notice that because I abcd depended only on e and DCI cross-ratios, the integrals 
labcd are manifestly dual-conformally invariant! And so (3.12) provides a manifestly 
dual-conformal representation of any one-loop ratio function in J\f—4 SYM. 

3.4 Theories with Triangle Contributions 

We should briefly mention that the regulator (3.2) can be applied to any planar 
theory, not just M = 4 SYM. (Planarity being a consequence of the way the momenta 
p a are associated with region-momentum coordinates x a ; a possible generalization to 
the non-planar case will be discussed shortly.) 

In a more general planar field theory, one-loop amplitudes may also require 
contributions from 'triangle' and 'bubble' integrals in addition to the scalar boxes. 

5 The factor of n in (3.9) should be thought of as '2n/2\ arising from the fact that each two-mass 
hard box diverges like \ log(e) 2 (one-mass boxes counted twice by symmetry), and the quad-cut 
coefficients f 1 and / separately generate the tree-amplitude for each of the n sets of adjacent legs. 
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These integrals manifestly break dual-conformal symmetry, but at least the triangles 
can be regulated in the same way as the boxes, (3.2). Indeed, regulated triangle 
integrals can be obtained from the 'DCI'-regulated scalar-box integrals without any 
additional work: one need only to send one of the points x a , . . . ,Xd to (space-like) 
'infinity', denoted Xoo. The correctness of this is easily seen from the geometry of the 
four-mass integral, (2.3). Thus, for example, the (no-longer 'DCT, but) e-regulated 
two-mass triangle integral can be found by simply sending a point of the three-mass 
integral to infinity: 

-I e a ,a + i,c, d = 14(1 - v) + \ logK) log(i;) + \ log(e) \og{v) + 0(e), (3.13) 
where we simply take X& — >%oc so that, 
'= lim /( a_1 ' 6 )( a ' &+1 )( c ' c? )l _ (a-l,6)(a,6+l) _ f (b,c)(a,d) \ _ (6, c) 
~~ d-*oo\(a-l,b+i)(a, c)(6, d)j (a- 1, 6+1) (a, c) ' ~~ d^oo\(a, c)(6, d)j (a,c)' 
Here, dual-conformal invariance is broken explicitly by the fact that d— >oo' picks out 

a preferred point — namely, x^ — in region-momentum space. The so-called 'bubble' 
and rational terms are unaffected by the infrared regulator. 

Our claim is that in any planar quantum field theory requiring triangle-contributions 
(which are absent for J\f—4), these contributions are correctly reproduced by simply 
augmenting the box expansion with triangle integrals obtained from l e ahcd as de- 
scribed above, with coefficients fixed by on-shell diagrams. The inclusion of bubble 
and rational terms would be unchanged from their usual form (see e.g. [35]). 

3.5 Integrand-Level Infrared Equations and Residue Theorems 

The well-understood factorization structure of infrared divergences of one-loop am- 
plitudes in gauge theory [36, 37] leads to constraints on integral coefficients in the 
framework of generalized unitarity, resulting in the so-called 'infrared (IR) equations' 
(see refs. [33,34,38] for a few applications). These equations can be understood in 
terms of unitarity cuts, and in the planar case lead to n(n-3)/2 relations among the 
integral coefficients. 

By considering generalized unitarity at the integrand level, however, the fac- 
torization structure of soft-collinear divergences leads to more powerful identities. 
Recall that the box coefficients are so-called 'quad-cuts' (co-dimension four residues) 
of the loop integrand; let us now consider 'triple-cuts' — co-dimension three residues. 
Because the loop integral is four-dimensional, a triple-cut still depends on one in- 
tegration variable; by applying Cauchy's residue theorem to the integral over this 
remaining variable, we find that the sum of all box coefficients sharing a triple-cut 

must vanish 6 . 

6 We should mention that there can be residues of poles at infinity — the parity-even combinations 
of such are simply the so-called triangle coefficients; the inclusion of these terms would have no 
substantive effect on the results of this section. 
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The richest of these residue theorems arise for triple-cuts involving at least one 
massless corner, as these separate into two distinct classes depending on whether the 
3-particle amplitude at the massless corner is A z or A^K Importantly, Cauchy's 
theorem applies to these two cases separately, leading to a pair of identities: 

a-1 d-l 

Ja,b,d,e J a,b,c,d ~ ia,b,d^n , 

e=d+l c=b+l 
a-1 d-l 

?2 _ _ A (k),0 



(6 = a+l) (3.14) 



fa,b,d,e fa,b,c,d ~ T a,b,d «Ai 

e=d+l c=6+l 

where r a) 6 j{ f = ±l if d = b+l or c? = a — 1, respectively, and vanishes otherwise. The 
sums appearing in the left-hand side of (3.14) are over all box coefficients sharing 
a particular (chiral) triple-cut, while the right-hand sides follow from the universal 
structure of soft-collinear divergences (and when r^O, (3.14) simply represents the 
famous 'BCF' formula for tree-amplitudes, [6,33]). 

It is easy to see that these constitute 2n(n — 3) linearly-independent equations. 
Averaging the two lines of (3.14) results in n(n — 3) relations among the (parity- 
even) box coefficients — which is twice as many as the "standard" IR-equations. This 
doubling of equations is a consequence of the fact that each collinear divergence of the 
integrated amplitude can arise from two distinct integration regions; and requiring 
factorization at the integrand level results in identities for each integration region 
separately. 

To better understand the nature of the vanishing terms appearing on the right- 
hand side of (3.14), let us consider a loop- integrand in the neighborhood of a soft- 
collinear divergence — where three consecutive propagators are simultaneously put 
on-shell. For the sake of concreteness, we may parameterize the loop integrand in 
such a region by writing t— > (L + , L~ ,£±) (see [39] for a similar discussion): 

X(£) ^ fc) '°/ ((L++ 1)L -_ p ±){L+L -_ tl){L+{L-+ 1) - t\) 1>{i) ■ (3 - 15) 
We can take a residue in L~ which puts the first propagator on-shell, resulting in 

^jw^t T ^ (3 ' 16) 

(for the regime in which L + , £j_ are both small). If we now write £\ =zz and drop 
the condition that z is the complex conjugate of z, we see that we can take two more 
residues — for instance L + and then z — to localize to the triple-cut, and then finally 
take a residue involving the pole in z. Notice that this results in a 'quad-cut' (a 
co-dimension four residue) which only involves three propagators! This is a simple 
one-loop example of the phenomenon of 'composite leading singularities' discussed 
in refs. [8,40], and is the physical origin of the terms on the right-hand side of (3.14). 
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It is not hard to show that the criterion for the convergence of an integral dis- 
cussed in section 3.1 is equivalent to the requirement that all composite leading 
singularities vanish (in general, composite leading singularities are in one to one 
correspondence with IR-divergent triangle topologies, which are in correspondence 
with these equations); this implies that for ratio functions in M = 4 SYM the right- 
hand sides of (3.14) must vanish. Moreover, it is possible to show that the difference 
between the 'DCI'-regularized boxes of (3.12), and more familiar — e.g. dimensionally- 
regulated — expressions for the box integrals is necessarily proportional to the the 
right-hand-sides of (3.14) — which provides an alternate proof of the equivalence be- 
tween the two regularization schema when applied to manifestly finite observables. 
We leave the details of this discussion as an exercise for the interested reader. 

In summary, the residue theorems (3.14) encode the physical fact that IR- 
factorization occurs at the integrand level, and constitute n(n-3) independent parity- 
even relations among the box coefficients — twice as many as the IR-equations arising 
after integration. These relations are general — that is, not specific to A/" = 4 SYM— 
provided that triangle coefficients are included where they are required. 

3.6 Applications to Non-Planar Theories 

Although we have largely focused our attention on planar theories, it is worth em- 
phasizing that the increased number of IR equations upon considering integrands 
instead of integrals is completely unrelated to planarity: the same result holds for 
non-planar theories as well — including gravity. This motivates a generalization of 
the integrand-level regulator (3.2) to the non-planar case, at least at one-loop order. 

As a simple illustration, consider the 4-particle amplitude in M = 8 supergravity. 
By the no-triangle property, [41] , it is a combination of at most three box functions. 
It is not hard to see that the absence of collinear divergences — the vanishing of the 
right-hand-sides of the analogs of (3.14) — uniquely fixes their relative coefficients, 
leading to the following form for the amplitude: 

^(o),i;Af-8 0CU F(p 1 ,p 2 ,p 3 ,p A ) + sF(pi,p3,p2,p 4 ) +tF(p 1 ,p 3 ,p 4 ,p 2 ), (3.17) 

where s, t, u are the usual Mandelstam invariants and the box integrals are labelled by 
the momenta coming into the vertices. It is not difficult to check that this matches the 
correct expression. As a less trivial example, consider the 6-graviton amplitude (with 
arbitrary heliciticies). In principle, this amplitude could involve any combination of 
195 boxes; but the integrand-level IR-equations of the preceding subsection shows 
that the amplitude can be expressed in terms of at most 120 combinations. That 
is, we find 75 linearly-independent constraints of the form (3.14), as opposed to 
the mere 25 constraints arising from the previously known IR-equations. It would 
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be interesting to explore if these relations could be used to obtain more compact 
analytic forms for one-loop graviton amplitudes. 



4. The Chiral Box Expansion for One-Loop Integrands 

As we have seen, while the familiar box expansion reproduces all one-loop amplitudes 
post-integration, it does not match the full structure of the actual loop integrand. 
These can be easily understood in the case of MHV loop amplitudes, where the 
actual MHV loop integrand only has support on two-mass easy boxes involving l\— 
with vanishing support on all quad-cuts involving £ 2 . However, because the scalar 
box integrals are parity- even, their integrands always have unit-magnitude residues 
for both quad-cuts. This is easy to understand, as scattering amplitude integrands 
are generally chiral, while the scalar boxes are manifestly non-chiral. 

In this section, we describe a slight modification to the scalar-box integrands 
given above which leads to a fully-chiral generalization of the box expansion, allow- 
ing us to represent all one-loop integrands of A/" = 4 SYM. In fact, such a modification 
was discovered for MHV and NMHV one-loop integrands in ref . [16] , but the general- 
ization to more complicated amplitudes was unclear. Here, by revisiting the special 
case of MHV, we will find that the underlying structure naturally generalizes to all 
N fc MHV one-loop integrands, A^' 1 . 

All of the essential structure needed is already present in the case of the MHV 
one-loop integrands, An^' 1 . In ref. [16], the one-loop MHV integrand was written: 



E (l (q-l a a+l)n(c-l c c+l))(Xca) (4V 
(£a-la)(£aa+l)(£c-lc)(£cc+l)(£X)' ' ; 

a<c<a — 9 9 — a<c<a 



Here, X is an arbitrary, auxiliary line in momentum-twistor space (of course, the 
integrand is ultimately, algebraically independent of X). We will mostly take this 
formula for granted here; but let us see what role is played by the auxiliary line X, 
and how we may generalize this to reproduce any one-loop integrand. 

Because a pentagon integral has five propagators, it has 2x(J =10 fourth-degree 
residues — from the two ways of cutting any four of the five propagators. But because 
one of its propagators involves this auxiliary line X, only two such residues are physi- 
cally meaningful: those which cut the lines {(a— la), (aa+1), (c-lc), (cc+1)}. These 
are obviously 'two-mass easy' quad-cuts; but because the numerator of the integrand 
in (4.1) is proportional to {£ (a-laa+l)f](c-l cc+l)), the pentagon's residue on £ 2 
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vanishes, while the residue on i\ is unity (see Table 3). But this is perfect: the 
one-loop MHV integrand only has support on £\\ 

Because of this, we choose to view the pentagon contributions to (4.1) as some- 
thing like a 'chiralized' version of the scalar two-mass easy box: 



T 1 

-^a,a+l,c,c+l 




(£ (a-1 a a+l)D(c-l c c+1)) (X c a) 
{£ a-la) (£ aa+1) {£ c-lc) (£ cc+1) {£ X) ' 



(4.2) 



We are motivated to draw this as a box because it has precisely one physically- 
meaningful quad-cut (in this case £±), upon which it has residue -1. Although not 
relevant for MHV integrands, we could similarly define the parity-conjugate version, 



X 2 




(£ca) (X (a-1 oa+l)D(c-l cc+1)) 
t a-la) (£ aa+1) {£ c-lc) (£ cc+1) (£ X) ' 



(4.3) 



which has residue of -1 on the quad-cut £ 2 , and vanishing residue on t x . 

Although it may seem like we're nearly done, we must step back to observe 
that not all the terms in (4-1) are pentagons! This is indeed a good thing, because 
as described in ref. [16], all such chiral pentagons are convergent, while the actual 
MHV one- loop amplitude is of course divergent! The easy-to-overlook, non-pentagon 
contributions to the one-loop MHV integrand of (4.1), come from the boundary terms 
when c = a+l: 



x; 



div 




a+\ 



(a-1 a a+l)n(a. a+1 a +2)) (X a+1 a) 
(£ a-la) (£ aa+1) {I aa+1) (£ a+la+2) (£ X) ' 

(£ aa+1) (a-laa+la+2) (X a+1 a) 
(£ a-la) (£ aa+l)(£ aa+1) (£ a+la+2) (£ X) ' 

(a-laa+la+2) (X a+1 a) _ 
(£ a-la) (£ aa+1) (£ a+la+2) (£X) 



(4.4) 



a+l 



a 
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We are motivated to draw this as a triangle, because it has only three non-X prop- 
agators. In fact, in the particular case where X is taken to be the point at infinity, 
X^ lv becomes precisely a scalar triangle. 

Therefore, although somewhat less concise than (4.1) (which is deceptively so), 
we can write the MHV one-loop integrand in the somewhat more suggestive form, 

^ 0),1 = E (5, a+ w + i) + E x " iv ' ( 4 - 5 ) 

l<a<c<n a 

where the divergent part is expressed solely in terms of triangles. (Although we have 
not yet defined all chiral boxes Tl bcd , the two-mass easy boxes given in (4.2) are the 
only ones relevant for MHV (k = 0).) 

From this simple example, it should be clear that if we have 'chiral' versions of 
all the scalar boxes — ones with precisely one physical quad-cut with residue -1 on 
either l\ and I2 — then would have a 'chiral' version of the box expansion: 



A(k),l _ V f f l 7-I , f2 j2 \ , A(k),0 V^jdiv 

•^r! — / , \J a,b,c,d-^a,b,c,d ' J a,b,c,d-^a,b,c,d J ' ^Si / . a • 



l<a<b<c<d 



(4.6) 



Notice that this expression will be valid before integration — unlike the more familiar 
scalar box expansion described in section 2. To see that this formula must be right, 
first observe that the chiral boxes with coefficients are specifically engineered to have 
precisely the same residues on all physical quad-cuts as the actual amplitude's inte- 
grand. However, as mentioned above, every chiral box integral I^l cd = jd 4: iX^ bcd 
is convergent — and so the chiral boxes alone cannot fully represent the amplitude. 
The remaining contributions must therefore encode the divergence of the one-loop 
amplitude, which from (4.5) together with integrand-level factorization is simply the 
sum of the divergent "triangles" 2^ lv . 

Thus, by construction, (4.6) will have the correct leading singularities on all 
physical quad-cuts (those that do not involve X) and have the correct infrared di- 
vergences. In order for (4.6) to reproduce the integrand, it obviously must be X- 
independent; and so it must have vanishing support on all quad-cuts involving X. 
The preceding arguments only fix the infrared-divergent part but not any potential 
four-mass integrals involving X. Thus, in order to uniquely fix X\ b c d and X^ bcd , we 
must also require that they are parity-odd on all four-mass quad-cuts involving the 
auxiliary line X — that is, they must have equal residues (in both sign and magnitude) 
on all 'spurious' quad-cuts (as any parity-even integrand must have opposite residues 
on parity-conjugate quad-cuts, this will ensure that no four-mass integrals involving 
X will survive integration over a parity-invariant contour). 

Given chiral box integrands i\ hcd and X a 2 bcd which satisfy the conditions de- 
scribed above — being infrared convergent, having only one nonvanishing physical 
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quad-cut and being parity-odd on spurious quad-cuts — it is not hard to prove that 
(4.6) matches all physical residues of the actual loop integrand, that it is ultimately 
free of any quad-cuts involving X, and that it is moreover algebraically independent 
of X. And so, (4.6) must give the correct integrand for the amplitude. (Using the 
Mathematica package documented in Appendix C, it is easy to verify that (4.6) di- 
rectly matches the one-loop integrands obtained using the BCF W recursion relations 
(which is described in Appendix B).) 

In terms of the chiral boxes, the ratio function becomes, 




(4.7) 



Notice that the divergent integrands X^ lv are manifestly canceled in the ratio function, 
leading to an expression involving only manifestly convergent chiral boxes. This 
is remarkable, as it provides an analytic form of any one-loop ratio function for 
which no regularization is needed! However, although the complete integrand is of 
course independent of X, each chiral box individually depends on X. Nevertheless, 
although it is generally difficult (algebraically) to prove the X-independence of the 
integrated expressions (this can be seen to amount to the integrand-level IR equations 
(3.14)), this remarkable fact can easily be verified numerically — for example, using 
the Mathematica package 'loop_amplitudes' which is made available with this 
note on the arXiv, and documented in Appendix C. 
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Table 4: Chiral One-Loop 'Box' Integrands and Integrals 




-\{Aa Cc)(BbDd){t{X))A 



-\{AaCc)(BbDd)(£ (X))A 

-ie«"^(Lon(^o(i*)n(i^2)) 



{la, L 2 , L 3 , L 4 } = {(Aa), (Bb), (Cc), (Dd)}, A = ^{\- U -vf-Auv 
u = u[ab;cd], v = u[bc;da], a = h(l—u + v+A), ft = |(l+u— v+A) 



1,2 . 



-Li 2 (a) -Li 2 (/3) + Li 2 (l) - \ log(w) log(v) +log(a) log(/3) 





a=B l((£(Aab)(M.Cc(DdB)r\(X))) 
-(£(Aab)f](Dd(CcB)f](X))) 



^(eB(Dd)n((Cc)f](Aab)(X))) 
-(iB(Cc)n((Dd)f](Aab)(X)))) 



I 1 ' 2 = Li 2 (l - ii[cfe;aX]) - Li 2 (l - u[bc;da}) 

+Li 2 (l -u[da;bX}) - \ log(u[cd;a X}) log(u[d a;b X}) 
— \ log(u[cb;aX]) log(w[d c;b X]) 




c C 




a=B l((eb(Aa)(M.Cc(DdB)r\(X))) 

-{£b(Aa)f](Dd(CcB)f](X)))) 



^(£B(Cc)n(Aa(Ddb)f](X))) 
-(£B(Cc)f](Dd(Aah)n(X)))) 



C=b 



I 1 ' 2 = Li 2 (l - u[dc;bX}) - Li 2 (l - u[db;aX}) 
+\ log(w[a b;cX]) log(u[db;aX]) 



1 



log(u[a c;dX]) log(u[d c;b X]) 




a=B 



C=b 



a=B 




{£(Aab)(\(C'cd)){XDB) 



(£DB){X{Aab)f](Ccd)) 



D=c 



I 1 ' 2 = Li 2 (l — u[db;aX]) — Li 2 (l — u[bc;d a] ) 
+Li 2 (l - u[ac;dX}) - Li 2 (l - u[cb;aX}) 
q +Li 2 (l — u[bc;dX]) — log(u[d b;a X]) \og(u[a c;d X]) 





a=B 



(£(Aab)r\(Ccd)){XDB) 



(£DB)(X{Aab)r\(Ccd)) 



D=c 



I 1 ' 2 = Li 2 (l) - Li 2 (l - u[db;aX}) - Li 2 (l - u[ac;dX}) 
r^-u —log(u[db;aX])log(u[ac;dX]) 



D=c 




(a- la a+1 a+2) (X a+1 a) 



(a- la a+1 a+2) (X a+1 a) 



a+1 



with 



rl,2. 



-Li 2 (l) — \ log(w[a-l a+l;a+2 X]) log(u[a+3 a+l;aX]) 



u[ab;cd] 



■|log(e) log(«[a-l a+l;a+2 X]u[a+3 a+l;aX]) — \ log(e) 

_ (a-la 6-16)(c-lc d-ld) 




(a-la c -lc)(b-lb d-ld) 



and ?i[a6;cX] 



a+1 

(a-la b-lb)(c-lc (X)) 



a 



(a-la c-lc)(b-lb (X))' 
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5. Conclusions 

In this note, we have revisited the familiar story of using generalized unitarity to 
reconstruct one-loop amplitudes, especially for the case of planar, M = 4 SYM. In 
order to make manifest all the known symmetries of the theory, we reconsidered the 
regularization of IR-divergences, and found a new, 'DCT-regularization scheme (3.2) 
which makes manifest — term-by-term — the dual-conformal invariance of all finite 
observables of A/" = 4 SYM at one-loop, including all N fc MHV ratio functions, (3.12). 

The existence of such a regularization scheme was motivated by considering the 
remarkable properties of one- loop amplitudes prior-to integration. Such considera- 
tions also led to integrand-level IR-equations, (3.14) — giving novel constraints among 
box coefficients, and having applications beyond the planar limit. And in section 4, 
we found that the familiar box expansion could be upgraded to a chiral box expan- 
sion for one-loop integrands, reproducing both the parity-odd and parity-even con- 
tributions to scattering amplitudes, and making the factorization of IR-divergences 
manifest. 

For the sake of completeness and reference, we have comprehensively described all 
the ingredients required to compute any one-loop amplitude in planar M = 4 SYM: we 
gave explicit formulae for all 'DCI'-regularized scalar box integrals in Table 1, and we 
gave expressions for all one- loop box coefficients in Table 3. Remarkably, these tables 
are incredibly redundant: all degenerate cases of both tables follow smoothly from 
the generic case — from the four-mass integral, (2.6), and the four-mass functions, 
(2.18) and (2.19), respectively. 

The notation used throughout this paper is reviewed in Appendix A. In Appendix B, 
we use the BCF W recursion relations described in [15] to explicitly represent all one- 
loop integrands in J\[ = 4 SYM. All the results described in this paper have been 
implemented in a MATHEMATICA package, 'loop_amplitudes'; instructions for ob- 
taining this package, and complete documentation of the functions made available 
by it are described in Appendix C. 

Natural extensions of this work include applying the 'DCI'-regulator to scalar 
integrals beyond one-loop, finding chiral representations of higher-loop integrands, 
and using the integrand-level IR-equations to find better representations of one-loop 
amplitudes beyond the planar limit. 
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A. Review of Momentum- Twistor Variables and Notation 



Momentum-twistor variables, introduced by Hodges in ref. [24], trivialize the two 
ubiquitous constraints imposed on the external kinematics for all scattering ampli- 
tudes: the on-shell condition, p\ = 0, and momentum conservation, ^2 a p a = 0. By 
this we mean that generic momentum twistor variables Z = [z\ • • • z n ) G G(4,n), 
with 2 a eC 4 , always correspond to a set of momentum-conserving, on-shell external 
momenta. This makes them especially convenient for use in scattering amplitudes. 

In section 2, we used region-momentum variables x a to encode the external four- 
momenta according to p a = x a+ i — x a . Momentum- twistors are so-called because they 
represent points in the twistor-space [25] of region-momentum x-space: 




The x-space polygon, whose definition depends on a choice for the cyclic ordering 
of the external legs, encodes the external momenta in a simple way. Each line in 
twistor space — spanned, say by the twistors (z a _i z a ) — corresponds to a point in x- 
space — in this case, the point x a . This is why, for example, integration over a point 
Xi in region-momentum space translates to integration over a line (£) = (£a^b) i n 
momentum-twistor space (see e.g. [11,15,16]). 

Given a set of momentum-twistors z a — viewed as columns of the (4 x n)-matrix 
Z — it is easy to construct the corresponding set of four-momenta. If we decompose 
each momentum-twistor z a according to 

Z a = (KfJ>a), (A.l) 

and define a (2xn)-matrix A = (Ai • • • A„) C (A -1 ) according to, 

r n n 5 a „ lb (aa+l) +5 ab (a+la-l) + 5 a+lb (a-la) 

X = fi-Q with Q ab = — — , (A.2) 

(a- 1 a) (a a+1) 

where (a b) = det{A a , X b }, then A-A = because Q-X = 0; as such, we may identify 
p a = A a A a , and these (on-shell) four-momenta will automatically conserve momentum. 

Conversely, given four-momenta written in terms of spinor-helicity variables ac- 
cording to p a = X a X a , momentum twistors z a can be constructed by joining each A a 
as in (A.l) with a corresponding /i a given by 

(ba) if 1 <b<a 



l_i = Q-X with Q 



ab - 







otherwise 



(A.3) 
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Supersymmetry is encoded by dressing each momentum-twistor z a with an anti- 
commuting four-vector r] a — collected into a (4xn)-matrix 77 acted upon by the 577(4) 
i?-symmetry of Af = 4 SYM. If we similarly define rj =rj ■ Q, then the kinematical 
data specified by {A, \,rj} will automatically be swpennomentum-conserving. 

Dual-conformal transformations in region-momentum space translate to mere 
5X(4)-transformations in momentum-twistor space; hence, dual-conformal invari- 
ants are written in terms of simple determinants: (abed) = det{z a , z^, z c , Zd}- The 
simplest dual-stzpenconformal invariant, however, involves five momentum-twistors, 
and is given by the familiar 5-bracket (sometimes called an '72-invariant'), [11, 12]: 
[abode] ~ ^^(Vaibcde) + Vb(cdea) + r) c (deab) + r] d (eabc) + r] e (abcd)) 

(abcd)(bcde)(cdea)(deab){eabc) 

All one-loop leading singularities except the four-masses 7 can be written directly 
as products of 5-brackets — as evidenced by Table 3 and the fact that BCFW recursion 
(see Appendix B) directly gives tree-amplitudes in terms of products of 5-brackets. 
These often involve geometrically-defined, auxiliary points in momentum-twistor 
space such as u (ab)f](cde)" which represents "span{z , zj,} f] span{z c , Zd, z e }" ■ All 
such objects are trivially found via Cramer's rule, which represents the unique iden- 
tity satisfied by any five generic four-vectors: 

z a (bcde) + Zb(cdea) + z c (deab) + Zd(eabc) + z e {abcd) = 0; (A. 5) 

from this, it is easy to see that 

(ab) P| (cde) = z a (bcde) + Zb(cdea) = — (z c (deab) + Zd{eabc) + z e (abcd)).(A.Q) 

A similar, geometrically-defined object which appears in the BCFW recursion of 
one-loop amplitudes (see Appendix B) is "(a be) f] (de/)", which simply represents 
"span{2 a , Zb, z c } H span{zd, z e , Zf}" — which, in this (projective) line in 

twistor space; for the sake of concreteness, we can always write this explicitly as: 

(abc)f](def) = (ab)(cdef) + (bc)(adef) + (ca)(bdef). (A.7) 

Finally, we should recall that the Jacobian arising form the change of variables 

from momentum-space to momentum-twistor space is the full Parke- Taylor MHV 

super-amplitude [42], , , ~ 

* 2X4 (A-^ 2X2 (A-A) , A8 . 



(12><23>-..<nl> ' 

which explains why (when written in momentum-twistors) all MHV amplitudes are 
simply An '° = 1, ensuring the dual-conformal symmetry of all amplitudes in planar 
Af = 4 SYM, [12]; and so, throughout this paper, An^' 1 should be understood as 
the color-ordered, single-trace contribution to the £-\oop integrand for the n-point 
N fe MHV scattering amplitude divided by (A. 8), and in units of (^ 2 iV c >/(167r 2 ))^. 



7 Rccall that the four-mass leading singularities require a prefactor of ipi 2 ( see equation (2.18)). 
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B. The BCFW Representation of One-Loop Integrands 



As described in ref. [15] (see also [11]), all £-\oop integrands for scattering amplitudes 
in planar M = 4 SYM can be found by the BCFW recursion relations. In terms of 
on-shell diagrams, the BCFW recursion relations correspond to: 




Being more explicit about the ranges for the terms involved, the recursion becomes, 

= ]T A^ <g) 4?*>*« + FL (A^ 1 ) . (B.2) 

n=n L +n R -2 BCFW 
k=k L +k R + l 

i=e L +i R 

Working in momentum-twistor space, the BCFW bridge operation corresponding 
to the shift z„-fz„ + az n _i, for n R >3, is given by: 

•A£? )4l (g) ^nnit = ^=a L ^ ■ • • , a-1, a )[1 a-lan-ln]Ai h ^(a, a,..., n-1, n), 

BCFW 

where a = (a a-1) f] (n-1 n 1) and n = (n n-1) f] (a-1 a 1); when rif> = 3 and riL — n—1, 
the bridge simply results in, 

^f®-^' - 4^(1,-., n-1). 

BCFW 

And so, the 'bridge' terms of (B.2) are fairly straightforward to compute in momentum- 
twistor space, and the operations involved are the same regardless of the loop-levels 
£l and £r of the amplitudes being bridged. (Of course, at tree-level, only the bridge 
terms contribute to the recursion; and so the discussion so far suffices to recursively 
compute all tree-amplitudes, An^' , in A^=4 SYM.) 

More interesting are the "forward-limit" contributions, FL(^4^ 2 1 ^' £_1 ). It is 
easy to see that (B.2) gives rise to i levels of nested forward-limits. As described 
in [11], determining which terms from the lower-loop amplitude are non- vanishing 
in the forward limit is generally difficult (even the number of terms which survive 
becomes scheme-dependent beyond one-loop). However, for one-loop amplitudes, 
we only need the forward-limits of trees; and as described in ref. [11], if the tree- 
amplitudes are obtained by BCFW , deforming the legs which are to be identified in 



- 29 - 



the forward- limit, then the terms which vanish are precisely those involving three- 
particle amplitudes on either side of the bridge. Therefore, the only non-vanishing 
contributions are: 



FL 



^ ) '°(e B ,l,...,n,£ A )) 

FL (4 fc L ='a°+i(^ l,...,a-l,a)[£ A £ B a-la n]A^°(a, a, . . . , n, t A ) 



n L ,n R >4 
k=k L +k R 



= A nL=aUl !> • • • > a-l,a)K[(l a-1 a); (1 n-1 n)]A%*>> Q (a, a,..., n-1, n, 

n L ,n R >4 
k=k L +k R 

where a = (a a-1) f] (£ a £b 1), n = (wn-1) H (£4^,8 1), and £= (£a^b) fl (^-1 n 1), and 
the "kermit" K[(abc); (def)} — written in terms of the line £=(£ a £b) — is given by: 

Kl{abc);{def)] S -*t {eab){ebc){eca){ede){fefWfd y (B.3) 
Putting everything together, the one-loop integrand for any amplitude is: 



4^(1,..., n) 
= 4^(1,..., n-1) 
+ ^ Ag L l{ L (l,...,o-l,o)[la-la7i-ln]^^(a,a,...,7i-l,n) 

n=nL+"B-2 

fc=fc L +A: iJ +l a= ( a a-1) f| 0-1 «-!), ™ = (nn-l) f| (a-1 a 1) 

71=71 £+71^ — 4 ^ 

^,n fl >4 a = (aa-l)n(^sl), n = (n n-1) fl (^s 1), t = {t A l B ) fl (n-lnl) 

Here, the first two lines represent 'bridge' contributions — identical in form to the tree- 
level recursion — while the last line represents the forward-limits. Notice the striking 
similarity of the roles of the 5-bracket [1 a-1 a n-1 n] in the bridge-terms and the 
'kermit' if[(la-la); (ln-ln)] in the forward-limit terms. Indeed, the forward-limit 
terms can be understood as unitarity-cuts which are "bridged" by the kermit. 

This analysis can of course be continued to higher loop-orders by repeatedly 
substituting the structure above into the forward-limit contributions appearing in 
(B.2); this results in higher-loop "kermits" which can similarly be understood as 
'bridging' amplitudes across a unitarity cut. However, as our present work requires 
only one-loop integrands (and as the the complexity involved in the higher-loop 
'kermits' is considerable), we will leave a more general discussion to future work. 
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C. MATHEMATICA Implementation of Results 



In order to make the results described in this paper most useful to researchers, we 
have prepared a Mathematica package called 'loop_amplitudes' which imple- 
ments all of our results. In addition to providing fast numerical evaluation of loop 
amplitudes and ratio functions, the loop_amplitudes package also serves as a reli- 
able reference for the many results tabulated above (as any transcription error would 
obstruct numerical consistency checks). 

The package together with a notebook illustrating much of its functionality are 
included with the submission files for this paper on the arXiv, which can be obtained 
as follows. From the abstract page for this paper on the arXiv, look for the "down- 
load" options in the upper-right corner of the page, follow the link to "other formats" 
(below the option for "PDF"), and download the "source" files for the submission. 
The source will contain 8 the primary package (loop_amplitudes .m), together with 
a notebook (loop_amplitudes_demo .nb) which includes many detailed examples of 
the package's functionality. 

Upon obtaining the source files, one should open and evaluate the Mathemat- 
ica notebook 'loop_amplitudes_demo .nb'; in addition to walking the user through 
many example computations, this notebook will copy the file loop_amplitudes .m 
to the user's ApplicationDirectory [] ; this will make the package available to run 
in any future notebook via the simple command "<<loop_amplitudes .m": 

in[i]:= <<loop_amplitudes .m 



Out[l] : = 




Occasionally, the "source" file downloaded from the arXiv is saved without any extension; this 
can be ameliorated by manually appending ".tar.gz" to the name of the downloaded file. 
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C.l Glossary of Functions Defined By the Mathematica Package 
Abstract Representations of Box Expansions 

• boxCoef f icient[n_,k_,cut_:0 ] [abed ]: returns the scalar box coefficient /abed 

for the amplitude An 1 if cut is (its default); alternatively, it will give /^ bcd or 
/abed if cu ^ i s se t to 1 or 2, respectively. 

• chiralBoxExpansion[ii_,.k_]: returns an abstract (symbolic) list of contributions 

to the one-loop amplitude An' 1 in terms of on-shell diagrams and chiral boxes. 

• chiralBoxRatioExpansion[ii_,k_]: returns an abstract (symbolic) list of contri- 

butions to the one-loop ratio function TZn 1 =A^' 1 -A^' An^' 1 written in terms 
of on-shell diagrams and chiral boxes. 

• onlyCyclicSeeds[expressioii_] : returns the parts of a symbolic expression for a 

loop integrand or integral which generate the entire expression upon summing 
over cyclic classes. 

• scalarBoxExpansion[n_,.k_]: returns an abstract (symbolic) list of contributions 

to the one-loop amplitude A^' 1 in terms of on-shell diagrams and scalar boxes. 

• scalarBoxRatioExpansion[n_,if_]: returns an abstract (symbolic) list of contri- 

butions to the one-loop ratio function T^n 1 =A^' 1 -A^'°A^' 1 written in terms 
of on-shell diagrams and scalar boxes. 

(Internal) Symbols for Objects Used by loop_amplitudes 

• chiralBox[cut_][iegList_]:an undefined object that is used internally by the 

package to represent the chiral-box J cut — for cut either 1 or 2 — involving the 
corners whose external legs are specified by legList. To be clear, given legList 
expressed in the form legList = {{a, . . . , B}, . . . , {d, . . . , D}}, chiralBox would (ab- 
stractly) represent the chiral box (see Table 4). 

• onShellGraph[cut_][IegList_,kList_]:an undefined object that is used internally 

by the package to represent the on-shell diagram involving the corner amplitudes 
of type N k MHV with kekList, and involving the legs specified by legList. To be 
clear, given legList = {{a, . . . , B}, . . . , {d, . . . , D}} and kList= {k a , . . . , k d }, then 
onShellGraph would (abstractly) represent the on-shell diagram whose corner- 
amplitudes are given by, 

{ At )fi (a h2 , a, . . . , B, fo), Xl b) ' (/3i,2, b, . . . , C, 71,2), fG . 

^ ( 7li2 ,c,...,D,M, A^ ) '°(5 1)2 ,d,...,A,a 1 , 2 ) }, 

for the quad-cut solution £± or £ 2 depending on whether cut is 1 or 2, respectively. 
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• Rlabcde ]: represents an undefined object that is used internally by the package 

to represent the 5-bracket involving twistors given by the sequence abcde. 

• scalarBox[IegList_ ] : an undefined object that is used internally by the package to 

represent the scalar box 7 abcd whose corners involve the external legs specified 
by legList, given in the form: legList = {{a, . . . , B}, . . . , {d, . . . , D}}. 

Analytic Expressions for Objects Involved in Box Expansions 

• boxRs[IegList_] : returns the list {Ri,R 2 } of the 5-bracket prefactors — for {£i,£ 2 }, 

respectively — (obtained as an on-shell graph where all corners other than those 
involving A% ^ have been replaced by MHV amplitudes) which multiplies the 
corner amplitudes to result in a given one-loop leading singularity. These are 
listed in detail in Table 3. 

• chiralBoxIntegral[iegList_] : returns the (manifestly finite) integrated form of 

the chiral box whose corners involve the legs legList. (Recall that both chi- 
ral integrands {I 1 ,! 2 } — being parity conjugates of one another — integrate to 
J 1 ' 2 because loop integration is always performed on a parity- invariant contour.) 
These integrals are listed in Table 4. 

• chiralBoxIntegrands[IegList_]: returns a list of the of chiral integrands {X 1 ,! 2 } 

for the box whose corners involve the legs specified by legList; the numerators 
of these integrands are listed in Table 4. 

• quadCuts[IegList_]:for a box whose corners are given by the legs specified by 

legList, quadCuts returns {{«i, . . . , 5\}, {a 2 , . . . , #2}}, specifying the points along 
(Aa), . . . , (Dd) which lie on the quad-cuts t\ and £ 2 , respectively; see Table 2. 

• scalarBoxIntegral[iegList_,JceepEpsiioiisQ_.'True ]: returns the integrated form 

of the scalar box whose corners are specified by legList using the DCI regulator 
described in section 2.2 (see Table 1). If the optional second argument of keep- 
EpsilonsQ is set to False, then all the divergences — proportional to log(e) — will 
be excluded, returning only J fin . 

• scalarBoxIntegrand[IegList_] : returns the canonically-normalized scalar box in- 

tegrand for the box whose corners are specified by legList. 

Replacing Abstract Objects with Analytic Formulae 

• fromRform[n_ ][expression_]: converts any 5-brackets in expression (encoded ab- 

stractly as R[a, b,c,d,e]) into explicit superfunctions encoded by pairs {/, C}. 
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• integrand[expression_] : converts any abstract boxes (scalarBox or chiralBox) oc- 

curring in expression to explicit integrands involving the loop-momentum {€)— 
represented by the line (£a^b) in momentum-twistor space. 

• integrate[expression_] : converts any abstract boxes (scalarBox or chiralBox) oc- 

curring in expression to explicit, integrated expressions. 

• xFoTia[expression_ ]: converts any abstract on-shell graph objects (onShellGraph) 

occurring in expression to explicit formulae written in terms of 5-brackets. 

BCFW-Recursion of Tree- Amplitudes and Loop Integrands 

• bcfwLoopIntegrand[n_,.k_]: returns the 'superfunction' form of rAmp[n,k, 1]- 

that is, not written in terms of 5-brackets, but in terms of pairs {/, C} which 
encode a superfunction f x5 kx4 (C -n). bcfwLoopIntegrand[n,i] is equivalent to 
fromRform[n] [rAmp[n,i,l]] . 

• rAmp[n_,]c_,£_:0 ]: returns the tree {(=0) or one-loop (( = 1) amplitude A^'^ as 

obtained via BCFW recursion — written in terms of 5-brackets and 'kermits' as 
described in Appendix B. 

• treeAmp[n_,i_]: returns the 'superfunction' form of rAmp[n, k, 0] — that is, not 

written in terms of 5-brackets, but in terms of pairs {/, C} which encode a super- 
function /x<5 fcx4 (C -rj). treeAmp[n,ic] is equivalent to fromRform[n][rAmp[n,if]]. 

Box Expansions of Loop Amplitudes and Ratio Functions 

• chiralLoopAmp[n_,l_] : returns a list of superfunctions contributing to the one- 

loop amplitude JdHAn^' 1 obtained using chiral boxes. chiralLoopAmp[n,ic] is 
equivalent to analyticIntegral[chiralBoxExpansion[n,]c]]. 

• chiralLoopIntegrand[n_,i_]: returns a list of superfunctions contributing to the 

one- loop amplitude integrand A^' 1 obtained using chiral boxes. chiralLoopIn- 
tegrand[n,if] is equivalent to analyticIntegrand[chiralBoxExpansion[n,i]]. 

• chiralRatioIntegral[n_,]f_]: returns a list of superfunctions contributing to the 

one-loop ratio function Jd 4 £ IZ^' 1 obtained using chiral boxes. chiralRatioInte- 
gral[n,i] is equivalent to analyticIntegral[chiralBoxRatioExpansion[n,i:]]. 

• loopAmp[n_,if_] : returns a list of superfunctions contributing to the one-loop am- 

plitude JdHAn ' 1 obtained using scalar boxes. loopAmp[n,.k] is equivalent to 
analyticlntegral [scalarBoxExpansion [n,k] ] . 

• ratioIntegral[n_,if_] : returns a list of superfunctions contributing to the one- 

loop ratio function Jd 4 £TZ^' 1 obtained using scalar boxes. ratioIntegral[n,i] is 
equivalent to analyticIntegral[scalarBoxRatioExpansion[n,]f]]. 
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Kinematical Specification and Numerical Evaluation 

• evaluate[expression_]: numerically evaluates all superfunctions occurring in ex- 

pression for the kinematical data defined by the global, (n x 4) matrix Zs. If 
expression involves an auxiliary line (X) or a loop- variable (£) = (£a^b), these 
are taken to be given by the last four entries of the global matrix Zs — that is, 

Zs= {Zi, ... ,Z n , Z Xl , Zx 2 ,Z£ A , z iB }. 

• exampleTwistors[n_] :it is sometimes convenient to evaluate analytic expressions 

using explicit kinematical data; under such circumstances, there are some con- 
veniences afforded by using "well-chosen" kinematical data. 

Reasons for preferring one choice over another include: having all Lorentz 
invariants be integer-valued and relatively small; having all dual-conformal cross- 
ratios positive (so as to avoid branch-ambiguities when evaluating the polyloga- 
rithms that arise in scattering amplitudes at loop-level); and possibly to have all 
Lorentz-invariants be distinct (either to help reconstruct an analytic expression 
or to avoid 'accidental' cancelations). Of these, the following momentum-twistors 
meet the first two desires spectacularly: 



The function exampleTwistros [16] is evaluated when the loop_amplitudes 
package is first loaded, allowing amplitudes involving as many as 16 particles to 
be evaluated without specific initialization. 

• randomPositiveZs[u_] : picks random kinematical data for which all cross-ratios 

are positive (the data Zs is positive when viewed as a four-plane: Zse £?+(4, n)). 

• setupUsingSpinors[iambdaList_,iambdaBarList_]:sets up the global variables 

Ls and Lbs for A and A, respectively, and defines the global (nx4) matrix Zs for 
momentum-twistors for use in numerical evaluation. 

• setupUsingTwistors[twistorList_]:sets up the global (rax 4) matrix Zs encoding 

the momentum-twistor kinematical data, and defines the auxiliary variables Ls 
and Lbs for A and A, respectively. 

• showTwistors: returns a formatted table illustrating the kinematical data — as en- 

coded by the currently used ones for evaluation by evaluate [] . 



Zs 



/l I I I 

2 3 4 5 

3 6 10 15 






- 35 - 



• superComponent[comporiei!t ] \_superFunction_ ] : in the loop_amplitudes pack- 
age, a superFunction is always represented by a pair {/, C} — an ordinary function 
f(l, . . . , n) of the kinematical variables times a fermionic ^-function of the form, 

s M (c-v) = f[\ (av^vlr-viX (G-3) 

7=1 ^ai<---<a fc J 

where C is an (n x fc)-matrix of ordinary functions, and for each a — 1, . . . , n, T) a 

is a fermionic (anti-commuting) variable. To be clear, we consider each particle 

Grassmann coherent state [39] of the form, 
1 1 

| a) = \a) {} + r 1 I a \a) {I} + -r) I a r)l |a) {J)J} + ^vivtVa \ a ){i,J,K} + V 1 aVlvlv a l a >{i,2, 3 ,4} ! 
and if we use r a to denote the i?-charge of the a th particle according to, 



field 


helicity 


i?-charge(r a ) 


short-hand for r a 


!«){} 


+ 1 


{} 


P 


\a){i} 




{1} 


p/2(^ {4}) 


\ a ){i,J} 





{I, J} 




\ a ){I,J,K} 


1 

2 


{I, J,K} 


m/2(^{l,2,3}) 


«){1,2,3,4} 


-1 


{1,2,3,4} 


m 



then superComponent[r!,. . . ,r n ][superFunction] returns the component function 
of superFunction proportional to, n 

n n (°- 4 ) 

a=l l£r a 

— that is, the component-function involving the states |l) n • • • \n) Vn . 

General Purpose and Functions and iEsthetical Presentation 

• nice\_expression_] : formats expression to display 'nicely' by making replacements 

such as ab[x- • • y]i— >(x ■ • -y), a [1] i — >- Qfi, etc., by writing any level-zero matrices 
in MatrixForm, and by drawing figures to represent abstract objects given by 
objects such as onShellGraph or scalarBox. 

• niceTime[timeJnSeconds_]: converts a time measured in seconds timelnSeconds, 

to human-readable form. For example, 

in[i]:= niceTime[299 792 458] 

niceTime [3 . 1415926535] 
out[i]:= 9 years, 182 days 

3 seconds, 141 ms 

• random[objectList_ ] : returns a random element from (the first level of) objectList. 

• timed[expression_] : evaluates expression and prints a message regarding the time 

required for evaluation. 
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